The effects of meson mixing on dilepton spectra 
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The effect of scalar and vector meson mixing on the dilepton radiation from hot and dense 
hadronic matter is estimated in different isospin channels. In particular, we study the effect of 
a-Lu and p — ao mixing and calculate the corresponding rates. Effects are found to be significant 
compared to standard n-n and K-K annihilations. While the mixing in the isoscalar channel mostly 
gives a contribution in the invariant mass range between the two-pion threshold and the ui peak, 
the isovector channel mixing induces an additional peak just below that of the (j). Experimentally, 
the dilepton signals from p-ao mixing seem to be more tractable than those from a-u mixing. 
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bJO! I. INTRODUCTION 

• There exists physical processes which are forbidden in free space but can take place in matter. Those are basically 
' related to medium-induced symmetry breaking effects. The interaction respects the symmetry, but the ground state 

breaks it For instance, in matter, Lorentz symmetry is lost which leads to the mixing of different spin states 

even when the interaction Lagrangian respects all the required symmetry properties. A well-known example of this 
^ , is a-u! mixing in nuclear matter as was first shown by Chin Q or p-ao mixing as pointed out in Ref. jsj. In addition 
there might be mixing between different isospin states in asymmetric nuclear matter. The p-LU mixing provides one 
such example as suggested in Q. 
I Our goal here is to identify the effects of scalar and vector meson mixing on dilepton production rates, which 

■ could in turn be observed in high energy heavy ion experiments. It is well known that the electromagnetic radiation 
, provides a penetrating probe to study various in-medium properties of vector mesons. Stimulated by the idea that in 

nuclear matter the vector meson properties might be modified from their vacuum values as a precursor phenomenon 
, to chiral symmetry restoration, a great deal of effort has been directed to understand their properties in hot and/or 

• dense nuclear matter [^,0. Several experiments have measured, or are planning to measure, dilepton pair production 
in nucleus-nucleus collisions. Such measurements have been carried out by the Dilepton Spectrometer (DLS) at 
Lawrence Berkeley National Laboratory g, and by HELIOS § and CERES ^ at CERN. The NA50 collaboration 
also measures dileptons yields in the light vector meson mass range pi] ]. Two new initiatives to be mentioned in 
this context are PHENIX at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 

■ and HADES at the Gesellschaft fiir Schwerioncnforschung (GSI) [Q. Density-dependent effects can also be found 
^> , in experiments at the Thomas Jefferson National Accelerator Facility (TJNAF) |lj]. As the effects to be discussed 

5^ ' here are mostly density-driven, the latter two facilities represent the energy regime and physical conditions perhaps 
more relevant for the present study. However, the physics discussed here could also be highlighted in an eventual 
low-energy RHIC run. 

For the understanding of the data in these experiments it is important to uncover the various processes which might 
lead to the production of dileptons and to estimate their relative contribution to the total production rate. Some effort 
has already being put in obtaining some "standard" sources related to the decay channels of various light mesons and 
processes like tt — tt or K-K annihilation ||l5|^^17[| . But the opening- up of new channels related to scalar- vector mixing, 
which constitutes the material to be presented in this paper, has only been addressed recently. 

The mixing of mesons considered here includes both the isoscalar and isovector channels. To be more specific, we 
estimate the rate of dilepton production from a-w and p-ao mixing. Even though some attention was paid to the 
former in the context of heavy-ion collisions , the importance of the latter was shown only recently in a zero 

temperature estimate of the dilepton production cross section ||] . The studies of the effects of mixing made so far 
were limited to calculations of cross-sections. 

This paper reports on calculations of the thermal production rates of dileptons induced by scalar- vector mixing in 
the isoscalar and isovector channels. Those results are then compared against standard tt-tt, K-K annihilations, at 
different temperatures and chemical potentials. Note that processes involving mixing of different G-parity states like 
a-u or p-ao are allowed only in matter which is not invariant under charge conjugation. Therefore, tt -I- tt ^ e~^e~ 
through their coupling to nucleons via a and w or tt + 77 — > e'^e~ mediated by p-ao mixing, can take place, as 
will become clear later, only in matter with finite nucleon chemical potential. Hence a natural place to look for such 
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density-driven effects are experiments involving relatively low energy where one expects to have a relatively large 
chemical potential. 

The model used here is Walecka-like with the inclusion of the aq meson and its mixin g w ith the p meson. The role 
of oo in the context of nuclear matter calculation has also been discussed in Ref . js],^ j22| . To keep the calculation 
simple we take nuclear matter to be symmetric under p n but this can easily be generalized to the case where 
nuclear matter has different number of protons and neutrons. 

In section II we discuss the general outline of the formalism. Next in section III we discuss a-LO and p-OQ mixing. 
In section IV we present our results for various physical cases. We summarize and conclude in section IV. 



II. FORMALISM 



The nucleon-meson interaction Lagrangian can be written as 

Cint = gaixpalt' + gao'4"t'ao,aT''i> + gu>N Nlpl^ill^Uj'^ + (1) 

gpl^^^T"^ + ^V^^M-^^S'-lpi^ , (2) 

where ip, (j)„, (paoi P ^-nd lj correspond to nucleon, a, an , p and w fields, and rs are Pauli matrices. The values used 
for the coupling parameters are obtained from Ref. |23|| . From now on we use s and v to denote scalar and vector 
mesons, i.e. s = <t, Oq and v = uj, p. It is understood that — 0.0. 

The polarization vector through which the scalar meson couples to the vector meson via the n-n loop is given by 

n^(go, m - 2*5,5, / ^Tr [G{k)T^G{k + q)] , (3) 

where 2 is an isospin factor and the vertex for v-nn coupling is: 

= 7m - ^^M'^'/" ■ (4) 
G{k) is the in-mcdium nucleon propagator given by 

G{k)^GF{k) + GT{k) (5) 



with 



and 



n (U^ 5{ko-Ek) , 6{k^ + Ek) , 

GT{k) = + ^J_[-^^^^_^ + ^,(^,^,)^^ ] (7) 

where = -^/k^ + m*^ and v = p — is the effective chemical potential. In Eq. (H) , the subscripts F and T refer 
to the free and temperature-dependent part of the propagator. Gf is the free Dirac propagator with the nucleon mass 
replaced by its effective value while the second term (Gt) includes the Fermi distribution of the nucleon assuming 
thermal equilibrium at temperature T. The delta function guarantees that on-mass-shell levels are occupied by the 
appropriate distribution functions, to* denotes the effective nucleon mass evaluated at the mean field level [^. It's 
value is given by the equation: 

where nk{T) = mE'^-u) — i nk{T) — f,(E'\,.) — ^^re the nucleon and anti- nucleon thermal distributions. It should 
be mentioned that in asymmetric nuclear matter the mean field generated by the qq meson would modify the neutron 
and proton mass differently. In the present case this effect vanishes. 
In this notation the nucleon density is written as 
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PB = ^ j cPk[nu{T) - nk{T)] . (9) 

Now we have a set of two coupled equations, Eqs. and ^) and two unknowns: the effective nucleon mass m* and 
the effective chemical potential v. One looks for the self-consistent solutions of these two coupled equations in order 
to find m* and v. The results are shown in Fig. f for two different nuclear densities. At a fixed density the effective 
nucleon mass increases and then again drops with temperature (T). This is consistent with Rcf. p^ . It might be 
mentioned here that we do not include resonances in the present work. 




T (MeV) T (MeV) 



FIG. 1. The in-medium nucleon mass and chemical potentialas function of temperature for two baryonic densities. 

Even though Eq. (^) involves Green's functions having both vacuum and temperature dependent part, the vacuum 
piece for 11^, as expected vanishes on account of Lorentz symmetry. Hence only the thermal part contributes to 
the mixed-polarization function and it has a natural cut-off provided by the thermal distribution functions. So, 
interestingly, the diagrams involving mixing of different spin states in matter is free from ultra-violet divergences. It 
should also be noted that the argument of 11^ involves go and \q\ indicating Lorentz symmetry breaking in matter. 
In general, it could be a function of q but for an isotropic medium it can only depend on the magnitude of the 
three-momentum. 

Returning to Eq. (||) , after the evaluation of the trace and a little algebra, 11^ could be cast into a suggestive form: 
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This immediately leads to two conclusions. First, it respects the current conservation condition g'^IIp = = Iliyg". 
Secondly, there are only two components which survive after the integration over azimuthal angle. In fact this 
guarantees that it is only the longitudinal component of the p meson which couples to the scalar meson while the 
transverse mode remains intact. Furthermore, current conservation implies that out of the two non-zero components 
of n^, only one is independent. After performing the angular integration Hq takes the form: 



f 



(11) 



From Eq. (^_I|) it is clear that the mixing amplitude vanishes for chemical potential v = 0. This is expected, as 
such mixing between a and lu (or p-ao) opens up new channels which violate G-parity conservation. This can only 
happen when the ground state is not invariant under charge conjugation as is the case in nuclear matter with finite 
nucleon chemical potential. In fact, isospin eigenstate u, which (ignoring small mixing with p induced by the quark 
mass difference and electromagnetic effects ) in free space cannot decay into two pions, in a thermal bath through 
its coupling with nucleon-antinucleon (or hole) excitation can decay into a pion pair with intermediate a meson [Q. 
The same argument holds true for p which in matter can decay into vr — 77 |q| . 

Note that the mixing amplitude vanishes at jgl = . Physically this corresponds to the situation when there is 
no way for the vector particle to distinguish between longitudinal and transverse direction and therefore rotational 
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invariance forces them to be degenerate forbidding mixing of different spin states in this limit. On the other hand, 
it has also been observed that at very high momenta Ho falls off. This corresponds to the situation when collective 
effects, which basically related to the long wave length oscillations, die off with high momenta. This is a well-known 
phenomenon in many-body physics. 



A. Isoscalar channel : a-u mixing 



As has been mentioned before that a-u mixing has been studied previously in the context of dilepton production 
for TT-vr annihilation in the s-wave channel in Ref. Jl8| , p^ . The thermal production rates of the dilepton from this 
channel however not been estimated before. In order to calculate the rate of dilepton production, wc first calculate 
the cross-section for the process Tr + Tr^CT^tj— >e+ + e^. The Feynman diagram for the relevant process are shown 
in fig |. 




FIG. 2. The Feynman diagram for the process tt -I- tt — » e+ -I- e 



The interaction Lagrangian for the ct-tttt vertex is given by 



^ = ^ffo-TTii-TOTr^cr'^Tr ' (fn (12) 



The cross-section can easily be evaluated to be : 
where, XIq is given by Eq. (pT|), and 



{AP - mlY + mlTl ' ' {M^ ^ mlf + mlTl{M) 



v/(Af 2 - Ami) 



For Ld we used a constant width of 9.5 MeV and for the sigma meson we adopt the same mass and width as Ref. 
p9[ . Recent discussion on this issue may be found also in Ref. |^,^ . We also have investigated the effects of form 
factors by considering the same form as in Ref. |18|: those are normalized to one for on-shell particles. Effects found 
are small as far as the dilepton yields are concerned. 



B. Isovector channel : p-ao mixing 



Quite similar to u-lj mixing p-a^ mixing opens up a new channel for dilepton production viz. vr + 
dense nuclear matter through n-n excitations. The Feynman diagram of the process is shown in Fig. 




FIG. 3. The Feynman diagram for the process tt + 77 ^ e"*" -|- e 
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To describe the iraQT] vertex we use 

2 ? 

For later convenience we define gTra„ri ~ faoTrrii'iTi^^ — rn^)/m-„. 

We should mention here that there is an uncertainty involved with the coupling parameter faa-nr; as discussed in Refs. 
p7| , p8| . This arises from the fact that oq lies close to the opening of the KK channel leading to a cusp- like behavior 
in the resonant amplitude, therefore a naive Breit-Wigner form for the decay width is inadequate. Furthermore, as 
mentioned before, there is also uncertainty involved with the aoNN coupling which renders the precise extraction of 
ao — TT — rj coupling even more difficult We take /aor;ir=0.44 which gives rap_,^^(mao)=59 MeV, while the 

experimental vacuum width of oq is between 50 — 100 MeV |^^. This choice is thus a conservative one. 

As before, the cross-section involving mixing in the isovector channel can be written in terms of momentum and 
temperature dependent mixing amplitude: 

_ 2TTa^ gl^^r, Fp{M^)Ga,{AP) 2 



where 
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{AP-mir + mlTjiM) ' {AP-mir+mlTUM) 



and the decay widths of the ao and p mesons: 



f2 2 _ 2 , 2 (M^ 2Y 



2M 



As was the case for the uj, we choose the mass of the aq according to Ref. |2^ ]. 



III. DILEPTON PRODUCTION RATES 



Once we know the relevant cross-section for the dilepton channel, the thermal production rate of the lepton pairs 
(number of reactions per unit time, per unit volume -Rf2'^ — dN^+^-pairs/d^^) can be estimated in the independent 
particle approximation of kinetic theory ||l5[| : 



dAP J (27r)3-"' "J {2Trf' dAP 



/l(ki) / ^f,(k2)—^is,M')Vrel (14) 



where /i, /2 are the thermal distributions of the 1,2 species and Vrei — '2eIe^ with \{x, y, z) = a;^ + + 

^ — 2xy — 2xz — 2yz is the triangle function. This could be cast into the form: 



\{A'P,ml,m'^) 



m (15) 



m = f Yi ^^^§^s[s - ip, + p^m;^- {ei, e2, ap) (w) 

'' i=l,2 ' 

where fi{Ei) = ^ob.^^ are the distribution functions for tt and rj mesons assuming thermal equilibrium. A departure 
from the standard notation should be noticed here. Usually, $(s) represents only the phase space factor |^^. In our 
case the invariant mass cross-section depends also on energies through Hq, we include that in the definition of $(s). 

After performing the angular integrations we get 

1 r^-^ + 

'^(*) = 7T^/ dEih{E^) dE2f2{E2)al^' {El,E2,Ad^) (17) 

{^T^r Jmi Je- 

where E^ = 2^[-E'i(s — m^) ± piX^^'^{s,mi,m2)]- <i>(s) is then calculated numerically. 
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A. Rates with a-u) mixing 



First we concentrate on the isoscalar channel or in other words the effect of a-LO mixing on the dilepton production 
rate. As the mixing amphtude now involves the u> propagator we observe a sharp peak in the vicinity of the uj mass. 
Obviously this mixing amplitude depends on the density and the temperature. The effects are more pronounced in 
Fig. H where we have higher density. No such peak is observed for the a meson which is sufficiently broad. We also 
plot the rates for tt-tt and K-K annihilation for comparison. The left panel shows the results at p/po — 1-5 and the 
right panel shows the same for p/po = 2.5. 




0.4 0.6 0.8 1.0 1.2 0.4 0.6 0.8 1.0 1.2 

M (GeV) M (GeV) 



FIG. 4. Effect of a-oj mixing on the dilepton production rate at T = 50 MeV for two different densities as mentioned in the 
legend. 

Fig.^ shows the results for the isoscalar channel at a higher temperature. It is evident that at higher density the 
effect of mixing is also high. 




0.4 0.6 0.8 1.0 1.3 0.4 0.6 0.8 1.0 1.2 

M (GeV) M (GeV) 

FIG. 5. Same as Fig. | but for T=75 MeV 

It is to be noted that to calculate the dilepton rate we have used tree level propagator for lj and a. In Ref. |T^ ] 
it is argued that a-u> mixing broadens the lj decay width substantially in nuclear matter at high density. This is 
understandable as now oj thorough nucleon-antinucleon (hole) excitations can decay into pion pairs as opposed to 
vacuum because of G-parity violation in matter. Therefore in matter this decay channel enjoys a larger phase space 
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and becomes as broad as p as its vacuum mass is almost degenerate with that of p meson. But at the one loop level the 
diagonal elements of the mixed propagator would also reduce its mass (off-diagonal elements have little effect on the 
effective mass) . Now if one considers the broadening induced by the off-diagonal element of the dressed propagator of 
the Lo meson, then to be consistent one should also take the medium modified mass for the uj meson. It is well-known 
that in Walecka model uo meson mass also drops p9| ] which would again reduce the phase space. These two competing 
effects are not expected therefore to change the in-medium width of uj significantly. We, however, adopt a different 
point of view. Our focus here is to see the effect of mixing on the dilepton production and how it compares with tt-tt 
annihilation. It should be noted that the contribution from a-uj mixing does not change the resultant spectra much 
expect a small peak representing uj meson pole in the transition amplitude. However, we expect an uj peak in the 
dilepton spectra even in the absence of meson mixing. The situation, as we shall see, is different for the isovector 
mode. There the dominant contributions are distributed near the tail of the tt-tt annihilation spectra (p spectral 
function). On the other hand if p spectral function gets sufficiently broadened this effect might become more visible 
near the lo peak region where additional contributions would come from the mixing. 



B. Rates with p-ao mixing 

Next we discuss the effect of mixing in the isovector channel which involves p-a^ mixing in Fig.| (|) for T=50 
MeV (75 MeV). This induces an ag bump just below the (p peak, p-a^ channel seems to be more interesting as it 
induces additional contribution beyond the tail of p spectral function. One can easily observe from Fig |^ that even 
at density p/po ~ 1-5 the contribution of the mixing wins over the tt-tt annihilation rate in the vicinity of M = 1.0 
GeV. Naturally at higher density this goes up as evident from the right panel of Fig. 
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FIG. 6. Effect of mixing on dilepton production rate in the isovector channel involving p-ao mixing at T=50 MeV. 

Interestingly enough the dilepton yield induced by the mixing in the isovector channel at higher temperature also 
seems to be significant as evident from the figures below. 
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FIG. 7. Same as Fig. | but for T=75 MeV 

Note that, unlike w, the broadening induced by the mixing on p decay width has been found to be marginal. One 
can understand this physically. Mixing opens up a new channel for p to decay into tt-t]. rj being quite heavy only the 
tail part of the p can contribute. Therefore it does not enjoy much of phase space to decay into tt — 77. It does seem 
reasonable to conclude that p-ap mixing provides more promising kinematical window in the invariant mass spectra 
to study exclusive in-medium effects. 

C. Rates with mixing both in isoscalar and isovector channels 

For completeness we show in Figs. |^ and ^ the total dilepton yield with the combined effects of mixing in both the 
isoscalar and isovector channels. It is important to see that rates induced by mixing are significant in certain windows 
of density and temperature. 
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FIG. 8. Rates at T = 50 MeV with mixing in both isoscalar and isovector channels. 



8 




0.4 0.6 0.8 1.0 1.2 0.4 0.6 0.8 1.0 1.2 

M (GeV) M (GeV) 



FIG. 9. Rates at T = 75 MeV with mixing in both isoscalar and isovector channels. 

IV. SUMMARY AND CONCLUSION 

To conclude and summarise, we have estimated the dilepton production rates arising from scalar and vector meson 
mixing in different isospin channel involving nucleon -hole (or anti-nucleon) excitations as intermediate states. The 
rates have been compared against the tt-tt annihilation in the p-wave channel (p-channcl) and K-K annihilation. 
Higher yield of dileptons in the invariant mass region between the p and the cf) peaks are expected. The mixing in the 
isoscalar channel contributes mostly between the two-pion threshold and the oj peak but is small compared to the tt-tt 
annihilation contribution. On the other hand, the isovector channel seems to provide a better probe as it contributes 
near the tail region of the two pion annihilation contribution. It should be stressed here that broadening or dropping 
p meson mass would even favor our observation. In that case the tt-tt background would be pushed more towards 
lower invariant mass region bringing the ao peak into a clearer relief. Finally, it is important to realize that mixing 
effects only appear in the longitudinal part of the vector meson spectral function and therefore there should be no 
effect of the mixing on the photon spectra which involves the transverse part of the vector meson polarization tensor. 

The effects studied here are density driven and can therefore be observed in environments where one has high 
baryonic chemical potential and comparatively lower temperature. The most relevant experiments to observe this 
kind of medium-induced phenomenon are the ones to be performed at GSI using the HADES detector where the 
density might reach up to three times the normal nuclear matter density and temperature T=60-90 MeV Ref. 
p2t also indicates a similar range of densities and temperatures at SIS/GSI energies. Those observations have guided 
our choice of the temperatures and densities used in this study. 

The calculations performed here can, and will be improved upon. It is important to proceed cautiously. As a next 
step, it would be necessary to address this problem with the full machinery of finite temperature many-body theory. 
This means, among other things, the dressing of the propagators and the consistent inclusion of baryonic resonances. 
Also, the interaction zone in a heavy collision is but a transient medium and the dynamics will contribute to mask 
our signal to some extent. However, owing to the nature of the electromagnetic signal, radiation from the dense 
regions will emerge. In our approach those dynamics would impact on the height of the oq peak. Similarly, a realistic 
assessment of the initial states of the dilepton-producing channels should involve a realistic simulation. A transport 
treatment can therefore address the concerns raised above. Such investigations are under progress. 
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